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ABSTRACT 

CN ; Based on RXTE/ASM and EXOSAT/ME data we studied X-ray variability of persis- 

tent LMXBs in the ^ 10^® — 10~^ Hz frequency range, aiming to detect features in 
their power density spectra (PDS) associated with the viscous time scale of the accre- 
tion disk ivisc- As this is the longest intrinsic time scale of the disk, the power density 
of its M variations is expected to be independent on the frequency at /< 1/tvisc- In- 
deed, in the PDS of 11 sources out of 12 we found very low frequency break, below 
^SJ ■ which the spectra are nearly flat. At higher frequencies they approximately follow the 

Pi, cx ly^^-^ law. 

, The break frequency correlates very well with the binary orbital frequency in a 

■ broad range of binary periods Porb ~ 12 min — 33.5 days, in accord with theoretical 
I expectations for the viscous time scale of the disk. However, the value of /break//orb is 

r~| . at least by an order of magnitude larger than predicted by the standard disk theory. 

CIh' This suggests that a significant fraction of the accretion M occurs through the optically 

Q I thin and hot coronal flow with the aspect ratio oi H/ R>Q.l. The predicted parameters 

$-H ' of this flow, T ~ 10'' — 10^ K and rig ^ 10^^ — 10^^ cm~"^ are in qualitative agreement 

, with recent Chandra and XMM-Newton observations of complex absorption/emission 

■ features in the spectra of LMXBs with high inclination angle. 

We find a clear dichotomy in the value of tvisc/^orb between wide and compact 

systems, the compact systems having ^10 times shorter viscous time. The boundary 
lies at the mass ration q « 0.3, suggesting that this dichotomy is caused by the 
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H excitation of the 3:1 inner Lindblad resonance in low-g LMXBs. 
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1 INTRODUCTION 

X-ray binaries show X-ray flux variations in a broad range 
of time scales (Fig0. Their power density spectra often re- 
veal a number of periodic and quasi-periodic phenomena, 
corresponding to the orbital frequency of the binary system, 
the spin frequency of the neutron star and quasi-periodic 
oscillation, which nature is still poorly understood. In ad- 
dition, aperiodic variability is observed, giving rise to the 
broad band continuum component in power density spec- 
tra, extending from ~ 10 — 100 msec to the longest time 
scales accessible for monitoring instruments. 

The aperiodic variability is often characterized by the 
fractional rms values of ~ tens of per cent, indicating that 
variations of the mass accretion rate M in the broad range of 
time scales are present in the innermost region of the accre- 
tion flow, where X-ray emission is produced. From the point 
of view of characteristic time scales, the high frequency vari- 
ations could potentially be produced in the vicinity of the 



compact object. Longer time scales, on the contrary, exceed 
by many orders of magnitude the characteristic time scales 
in the region of the main energy release. Hence, the low 
frequency Af variations have to be generated in the outer 
parts of the accretion flow and to be propagated to the re- 
gion of the main energy release. Due to the diffusive nature 
of the standard Shakura-Sunyaev disk dShakura fc SunvaevI 
1973), it plays the role of the low-pass filter, at any given 
radius r suppressing M variations on the time scales shorter 
than the local viscous time tvisc{r). The viscous time is, 
on the other hand, the longest time scale of the accretion 
flow. Hence, it is plausible to suppose that any given ra- 
dius r contributes to M and X-ray flux v ariations predom- 
inantly at the frequency / ^ 1/tvisc iLvuba rskii' 19971: 
IChurazov et alJl200lD . lLvubarskiil 119971) demonstrated that 
in this picture Pi, oc v~" power density spectra can naturally 
appear with slope a ~ 1, in good agreement with observa- 
tions. 



2 M.Gilfanov & V.Arefiev 




log(Frequency), Hz 



Figure 1. The power spectrum of Cyg X-2 in broad frequency 
range. Based on non-simultaneous ASM, PCA and EXOSAT 
data. The soUd hne shows = 1 ■ 10~* u~^'^ power law. 



Owing to the finite size of the accretion disk, the longest 
time scale which can appear in the disk is restricted by the 
viscous time on its outer boundary tvisc{Rd)- Below this fre- 
quency the X-ray flux variations are uncorrelated, therefore 
the power density spectrum of the accretion disk should be- 
come flat at /^/visc ~ 1 /tvisc{Rd) If there are several 
components in the accretion flow each having a different 
viscous time scale (e.g. geometrically thin disk and diffuse 
corona above it), several breaks can appear on the power 
spectra at the frequencies, corresponding to the inverse vis- 
cous time scale of each component. Yet another source of 
M variations might be the variability of the mass trans- 
fer rate from the donor star. As these variations also have 
to be propagated through the accretion disk, only low fre- 
quency perturbations, with / < l/fvisc(i?d), will reach the 
region of the main energy release. This might give rise to 
an independent continuum component in the power density 
spectrum at / < l/tvisc(i?d). The shape of the final power 
density spectrum of X-ray binary will depend on the relative 
amplitudes of the M variations due to different components 
of the accretion flow and the donor star. As these are inde- 
pendent from each other, a distinct features should appear 
in the power density spectrum of X-ray binary at frequencies 

/ ~ l/tvisc(i?d). 

The viscous time s cale of the accretion disk is 
JShakura fc Sunvaevlll97tl 



The broad QPO like features which can appear near t~^^^ due to 
the mechanism suggested bv [Vikhliniii ct al.. 1.1994) are discussed 
in section [5. 31 



t. - 2 (Hd\-' 1 n) 
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where Hd is the disk half thickness and Q,K{Rd) is the Kep- 
lerian frequency at the outer disk boundary Rd- For a geo- 
metrically thin disk, H/ R <^ 1, the viscous time scale is sig- 
nificantly longer than the Keplerian time at the outer disk 
boundary. If the disk fills a large fraction of the Roche-lobe 
of the primary, as is the case for steady disks in persistent 
LMXBs, Rd ~ 0.7-0.9 Rl, (sectionlOl. the viscous time at 
the outer boundary tvisc{Rd) will exceed the orbital period 
of the binary system Porb as well. Therefore, one might ex- 
pect to find features associated with the viscous time of the 
accretion disk at frequencies below the orbital frequency of 
the binary system. Detection of such features for a sample 
of LMXBs would give an opportunity to probe the accretion 
disk parameters, such as viscosity parameter a and the disk 
thickness at its outer boundary. 

In this paper we study a sample of persistent low 
mass X-ray binaries aiming to detect such features in their 
power density spectra. Longevity of typical orbital periods 
in LMXBs requires long monitoring observations to achieve 
this goal. This has become possible thanks to the long op- 
erations of the All Sky Monitor aboard Rossi X-ray Timing 
Explorer, which provided long term light curves of hundreds 
of galactic sources covering period of ~ 7 — 8 years. Some 
sources with rather short orbital period, < few hours may 
be studied based on the EXOSAT data. 

The paper is structured as follows. In section |21 we de- 
scribe the selection criteria and sample of LMXBs. Power 
density spectra and their approximations are presented in 
the sectional The viscous time scale of the accretion disks 
as obtained from the power spectra analysis and constrains 
on parameters of the outer disk are discussed in the section 
|1] In section|S]we consider the vertical structure of the semi- 
thick accretion flow, derive parameters of the disk corona 
and compare these with other evidence of the existence of 
diffuse coronal flow in LMXBs. Our results are summarized 
in section|H] In the Appendix we describe the method used to 
calculate power density spectra from the ASM light curves. 



2 THE DATA AND SELECTION CRITERIA 
2.1 Data 

We used data of the ASM instrument of RXTE observatory 
(Swank, 1999), covering the period from 1996-2004 (MJD 
range from « 50100 — 53100). The ASM instrument operates 
in the 2-12 keV energy range and performs flux measure- 
ments for over ~ 300 X-ray sources from the ASM source 
catalog once per satellite orbit, i.e. every ~ 90 min. Each 
flux measurement (dwell) has duration of ~ 90 sec. Due to 
navigational constrains and appearance of very bright tran- 
sient sources, the ASM light curves for individual sources 
sometimes have gaps of duration of up to few months. The 
dwell-by-dwell light curves were retrieved from the public 
RXTE/ ASM archive at HEASARC. 

For sources with short orbital periods, Porb ^ 4 hours, 
we used the data of the medium energy (ME) detecto r of 
EXOSAT satellite (iTurner. Smith fc Zimmermannll98lh . It 
provided ~ several tens of ksec long light curves with typical 
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Table 1. The binary system parameters for LMXBs from the ASM and EXOSAT samples. 



Source 


Ml 


M2 


^orb 


a 


q 


i 


D 


Rof. 




Mq 


Mq 


hrs 


10^° cm 




dog 


kpc 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


ASM sample 


GRS 1915+105 


14.0 + 4.4 


0.81+0.53 


804 + 36 


747 + 78 


0.058 + 0.033 


66° + 2° 


12.1 + 0.8 


1,2 


GX 13+1 


1.4 


5 + 1 


592.8 


461 + 25 


3.6 + 0.9 




7 + 1 


3 


Cir X-1 


1.4 


4 + 1 


398.4 


334 + 21 


2.9 + 0.8 




5.5 


4 


Cyg X-2 


1.78 + 0.23 


0.60 + 0.13 


236.3 


180 + 7 


0.34 + 0.04 


62.5 + 4° 


7.2 + 1.1 


5 


GX 349+2 


1.4 


1.5 + 1 


22.5+0.1 


40.0 + 4.7 


1.07 + 0.73 




9.2 


6 


Sco X-1 


1.4 


0.42 


18.92 


30.5 + 1.1 


0.30 + 0.05 


38° 


2.8 + 0.3 


7 



EXOSAT sample 



EXO 0748-676 


1.4 


0.45 


3.824 


10.6 + 0.4 


0.32 + 0.06 


~ 75° 


7.6 


8 


4U 1636-536 


1.4 


0.36 


3.793 


10.3 + 0.4 


0.26 + 0.05 




5.9 


9 


4U 1323-619 


1.4 


0.26 


2.93 


8.53 + 0.35 


0.19 + 0.03 


<80° 


15 + 5 


9 


4U 1916-053 


1.4 


0.09 + 0.02 


0.83 


3.55 + 0.16 


0.064 + 0.017 


>70° 


10.0 


10 


IH 1627-673 


1.4 


0.08 + 0.01 


0.690 


3.13 + 0.14 


0.057 + 0.011 


<8° 


8.0 11,12 




4U 1820-303 


1.4 


0.07 + 0.01 


0.190 


1.32 + 0.06 


0.050 + 0.010 




8.0 


13 



(1) - the mass of the compact object, if the uncertainty is not given a = 0.2Mq is assumed; (2) — the mass of the secondary, if the 
uncertainty is not given, 10% relative error is assumed; (3) - orbital period; (4) — binary separation computed from the 3rd Kepler law; 
(5) - mass ratio q = M1/M2. For GRS 1915+105 and Cyg X-2 from refs (1) and (5) respectively, for other sources computed from the 
masses given in the columns (1) and (2) of the table. (6) - inclination; (7) — source distance; (8) — reference for binary system parameters: 
1 Harlaftis & Greiner 2004; ^ Fender et al. 1999; ^ Bandyopadhyay et al. 1999; Johnston et al. 1999; ^ Orosz and Kuulkers 1999; ^ 
Wachter & Margon 1996; Steeghs and Casares 2002; * Parmar et al. 1986; ^ an estimate based on the empirical Po^;,— mass relation 
from Patterson, 1984; 1° Nelson et al. 1986; " Levine et al. 1988; Chakrabarty 1998; " Rappaport et al. 1987; " the companion mass 
estimate is very uncertain. 



time resolution of ~ 1 sec in the 0.9-8.9 keV energy range. 
The EXOSAT data were also retrieved from HEASARC. 



2.2 Selection criteria 

We have selected low mass X-ray binaries, satisfying the 
following criteria: 

(i) Persistent sources, which light curves do not show ofT- 
states or X-ray Novae-like outbursts. This ensures that the 
accretion disk was in the same state during the analyzed 
period. 

(ii) The orbital period Port is known. Orbital p arameters 
of X-r ay sources were taken from the catalog of iLiu et alJ 
J200lfl . 

The Nyquist frequency in the power density spectra ob- 
tained from the ASM light curves, /Nyq ~ 10"* Hz, corre- 
spond to period of ~ 3 hours. In practice, due to specifics of 
the ASM light curves, the high frequency part of the power 
spectra can sometimes be distorted by the aliasing effect. To 
minimize the contribution of this effect we selected source 
with sufficiently long orbital periods, Porb ^ 10 hours. 

Although the aliasing is not an issue for EXOSAT light 
curves, their duration, typically limited by ~ 50 ksec, im- 
poses a different constrains on the Poi-b- As we are looking 
for features at the time scales larger than Porb, we required 
that the light curve covers at least several Porb- This leads 
to the constrain for the sources which can be studied using 
EXOSAT light curves. Port < 4 hours. 

(iii) Source brightness. We selected sources with the aver- 



age count rate exceeding > 2 cnts/s for ASM and > 5 cnts/s 
for EXOSAT. 

(iv) The noise level, calculated for ASM power spectra, 
although approximately correct, is not accurate enough, 
due to existence of un accounted systemat ic errors in the 
flux measurements (e.g. iGrimm et al ]|2P02). To avoid addi- 
tional uncertainties in the power spectra we considered only 
sources with sufficiently high signal at the orbital frequency, 
Psignai ^ lOPnoisc at /orb . For EXOSAT light curvcs the noise 
level is not an issue. 

2.3 The sample 

There are 6 sources from the ASM catalog and 6 sources 
observed by EXOSAT, satisfying the selection criteria of 
the section 17!^ The final ASM sample includes the follow- 
ing sources: GRS 1915+105, GX 13+1, Cir X-1, Cyg X- 
2, GX 349+2, Sco X-1 The EXOSAT sample includes: 
4U 1323-619, 4U 1636-536, 4U 1820-303, 4U 1916-053, 
IH 1627-673, Exo 0748-676 

The binary system parameters for selected sources are 
listed in Table Q and their X-ray light curves obtained by 
ASM are shown in Fig. ^ 



3 POWER DENSITY SPECTRA 

Power density spectra were computed in the 2-12 keV 
(ASM) and 0.9-8.9 keV (EXOSAT ME) energy range. The 
PDS of the sources from ASM sample were obtained using 



4 M.Gilfanov & V.Arefiev 



^100 

(/) 

> 80 

-+-> 

o 60 
^ 40 
20 

o 60 

< 50 
a 40 

o 

S 30 

m 

< 20 

30 



< \ < \ < I I I 

ScoX-1 - 

^ • • • 

cnts/1 = 



I ' \ ' r 
•V: • CirX-l' 



J \ I 



^^I G'AS 191 '5+1 05 



• "• • • •• 



' GX344+2' - 



Cyg X^2 



Gi 13+1 



- ' 4 ♦ = 



4U 18^0-30a 



20 



^ 10 



4U 16:^6-53& 



t I »«- 



Ihi 16'27-'673 
cnts xlO 



4U '19^6-653- 

I 



4 U'l 3^3-61 9: 



1998 2000 2002 
Year 



1998 2000 2002 
Year 



Exb 0748-676 



\ I \ \ \ \ I 

1998 2000 2002 
Year 



Figure 2. Long term X-ray light curves of selected sources obtained by ASM (2—12 keV). The light curves cover the period from 
1996-2004 



the method based on the autocorrelation function calcula- 
tion as described in the Appendix |X] The EXOSAT hght 
curves were analyzed with the aid of the powspec task from 
FTOOLS 5.1. In analyzing the EXOSAT data we averaged 
the power spectra obtained in (nearly) all individual obser- 
vations with adequate time resolution available in the public 
archive at HEASARC. 

All but one source in the EXOSAT sample are X-ray 
bursters. The presence of X-ray bursts in their light curves 
can results in appearance of an additional component in the 
PDS, having no relation to the power spectrum of M varia- 
tions due to the accretion disk. To avoid this contamination, 



we screened the EXOSAT hght curves to exclude the time 
intervals corresponding to X-ray bursts. In the case of EXO 
0748-676 we also excluded the time intervals corresponding 
to the X-ray eclipses. No attempt to screen out the X-ray 
dips has been done, due to ambiguity of their identification. 
The power spectra based on the ASM data are not sub- 
ject to a significant contamination due to X-ray bursts and 
dips/eclipses. 

The power spectra of selected sources are shown in 
Figsl^and^] For reference, we plot in each panel a power 
law spectrum Pj, = 1 • 10"'* i^"^ "^. Flattening at frequen- 
cies / > 10~" Hz can be seen in the PDS of GRS 1915-H05. 
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Figure 3. Power density spectra of LMXBs from tlie ASM sample. Tlie solid line in each panel shows = 1 ■ 10"'* u~^'^ power law. 
The vertical dashed line marks the orbital frequency. In the upper left panel the histogram marked with hollow circles without error bars 
shows the power spectrum of GRS1915+105 obtained from the entire ASM time series of duration 8 years. Note that the vertical scale 
is different in different panels. 



6 M.Gilfanov & V.Arefiev 




Figure 4. Power density spectra of LMXBs from the EXOSAT sample obtained from the ME data. The solid line in each panel shows 
Pv = 1 ■ 10~*i^~^'^ power law, the same as in Fig|2l The vertical dashed line marks the orbital frequency. For EXO0748-676 and 
1H1627-673 the power spectra obtained from the GSPC data are also shown (histograms marked with hollow circles). Note that the 
vertical scale is different in different panels. 
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Figure 5. Broad band power spectra of LMXBs based on ASM and EXOSAT data. The solid line in each panel shows = 1 • 10^^ v^^-^ 
power law, the same as in Figl^andlH The vertical dashed line marks the orbital frequency. For Cyg X-2 the power spectrum obtained 
from the RXTE/PCA data is also shown (the histograms marked with open circles). Note that the vertical scale is different in different 
panels. 
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This can potentially be related to the specifics of ASM light 
curves, namely, it can be caused by the aliasing effect leading 
to the leakage of higher frequency power below the Nyquist 
frequency. The assessment of the reality of this flattening is 
beyond the scope of the present paper. Also seen for some 
sources (most prominently for Cir X-1 and 4U1820-303) are 
the peaks due to the orbital modulation. In Fig|^we show 
combined power spectra based on the ASM and EXOSAT 
light curves. These spectra cover ~ 10"® -10"^ Hz frequency 
range. For several faint sources from the EXOSAT sample 
no meaningful power spectra can be obtained from the ASM 
light curves due to too low count rate. For GRS1915+105, 
obviously, no EXOSAT data exist. 

The power spectra of almost all sources show a clear low 
frequency break, below which they are nearly flat. For some 
sources the second very low frequency component is present 
at lower frequencies. The most clearly this component can be 
seen in the case of Cir X-1, Sco X-1 and 4U1916-536. Some 
sources also show broad QPO like features near or below the 
break frequency. At high frequencies, above the break, the 
power spectra follow the oc power law. Remarkably, 
the slope of the power law appears to be similar in all 12 
sources, a ~ 1.3. The same is true for the normalization. 

3.1 4U 1636-536 and 4U 1820-303 

There is at least one exception from the above behavior. No 
evidence for break was found in the PDS of 4U 1636-536. 
Although the power spectrum of this source might have sev- 
eral weak features, its overall behavior in an extremely broad 
frequency range from 1Q~* — 10~^ Hz is well represented by 
the power law with slope a ~ 1.3 and without any obvious 
break (Figl^. 

The second possible exception is 4U 1820-303. In addi- 
tion to clearly visible peak due to the orbital modulations, it 
appears to have a shoulder at / ~ 10^'^ — 10~^ Hz followed 
by a power law at lower frequencies. The presence of this 
shoulder is more obvious in the data of individual EXOSAT 
observations, shown in FigEl In the formal statistical sense 
the existence of the low frequency break is highly significant. 
For example, approximation of the Aug. 19-20, 1985 data 
in the 10~^ — 2 • 10~^ Hz frequency range gave the follow- 
ing values: = 12-7 (11 d.o.f.) and = 42.2 (12 d.o.f.) 
for the power law model with and without low frequency 
break respectively, resulting in Ax^ ~ 30 for one additional 
parameter. However, considering the overall shape of the 
power spectrum (e.g. FiglSJ, this source presents a less ob- 
vio us case than the ot h ers in our sample. 

fc Grindlavl (|2Q^ analyzing long term X-ray 
variability of 4U 1820-303 confirmed earlier suggestion that 
it is a hierarchical triple system, where ultra-close X-ray bi- 
nary is orbited by companion with orbital period about 1.1 
day. Due to influence of the third star the well-known 176 
day X-ray modulation of 4U 1820-303 is generated. Influ- 
ence of the third star can also give rise to very low frequency 
modulation of the mass transfer rate in this system, unre- 
lated to the M variations intrinsic to the accretion disk. This 
could, in principle, explain strong very low frequency power 
law component at / < 10^'^ Hz. In the following we asso- 
ciate the shoulder at log(/) ~ —2.5 with the low frequency 
break similar to the ones observed in other LMXBs in our 
sample. We emphasize that this interpretation is not unique 




10-3 10-2 

Frequency, ilz 



Figure 6. The power spectrum of 4U1820-303 obtained by EX- 
OSAT on April 16 (open circles) and August 19-20 (filled circles), 
1985. Tlie solid line sliows = 3 • IQ-^ ;y-i-3 power law. The 
vertical dashed line marks the orbital frequency. 

and 4U 1820-303, like 4U 1636-536, can show the behavior 
different from other sources in our sample. 

3.2 PDS approximation 

The power density spectra were approximated by a model 
of a power law with low frequency break: 

i Af-" f> /break 

P{f) ^1 _ ■' - ■' (2) 

I ^ /break / < /break 

An additional Lorentzian component was added to the 
model for the sources with strong orbital modulation and/or 
broad low frequency QPO peaks. The frequency range was 
chosen for each source individually, depending on the du- 
ration of the available time series and the presence of the 
additional noise component below /break- The adopted fre- 
quency ranges and best fit values of the model parameters 
are listed in Table 2, along with the count rates and X-ray 
luminosities of the sources. The large width of the QPO 
component present in the power spectra of several sources 
close to the low frequency end makes the determination of 
the break frequency somewhat ambiguous. For these sources 
(GRS1915-f 105, Cyg X-2 and 4U 1323-619) we list the val- 
ues of the break frequency obtained with and without the 
Lorentzian component in the model. In the following the av- 
erage of the frequencies of the two models are used. The sta- 
tistical uncertainty of the break frequency in these sources 
are obtained from combined confidence intervals of the mod- 
els with and without Lorentzian component. At least three 
harmonics of the orbital frequency are present in the power 
spectrum of 4U1323-619. The corresponding bins of the 
power spectrum were excluded from the fitting procedure. 
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Table 2. The parameters of power spectra approximation. 
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cnts / sec 


cnts / sec 
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Hz 






Hz 
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(1) 
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GRS 1915+105 


60.7 


— 


34 


4.03+J;«^ X 10- 


8 
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— 


5 


10 
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- 1 ■ 


10- 


-6 










6.50ti;g4 X 10- 


8 




2.81:2:8 X 10"* 


5 


10 


-9 


- 1 ■ 


10" 


-6 
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23.1 


— 


4.3 


1.13 + 0.25 X 10 


-7 


0.56 + 0.03 


— 


1 


10 


-8 


- 1 ■ 


10" 


-5 
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— 


8.6 
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- 1 • 
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3.66l°;47 X 10- 


7 
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1.87 + 0.09 X 10-7 


1 


10 


-8 


- 1 • 


IQ- 


-5 


GX 349+2 


51.6 




16.7 


2.78+°-g2 X 10- 


6 


0.66 + 0.06 




1 


10 


-8 


- 1 • 


10- 


-5 


Sco X-1 


917 




27.4 


e.u+lf^ X 10- 


6 


O-Sy-0.06 




5 


10 


-7 


- 5 • 


10" 


-5 


EXO 0748-676 


0.57 


44.0 


0.37 


1.03lo;i2 X 10- 


4 


1.25 + 0.01 




1 


10 


-5 


- 3 • 


10" 


-2 


4U1636-536 


14.7 


153 








1.29 + 0.01 




1 


10 


-5 


- 3 • 


10- 


-2 


4U1323-619 


0.60 


4.4 


0.14 


2.01 + 0.22 X 10 


-4 






1 


10 


-5 


- 3 • 


10" 


-2 










2.10 + 0.60 X 10 


-4 
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(1), (2) — ASM and EXOSAT count rates; (3) — X-ray luminosity calculated using respective count rates and source distances listed in 
Table [Tl (4) - break frequency; (5) — power law slope; (6) - the centroid frequency of the Lorentzian component; (7) - frequency range 
for PDS fits. 



4 VISCOUS TIME SCALE IN LMXBS 

For the standard accretion disk (Shakura & Sunyaev, 1973) 
the viscous time scale at the outer edge of the disk Rd can 
be estimated as: 

where Rd is the disk outer radius, fix = GMi / R\ is the 
Keplerian frequency, Hd is the disk half-thickness at the 
outer edge. Mi is the mass of the primary and a is the 
dimensionless viscosity parameter. 

Combining the eq.J^ with the third Kepler law: 

Porb = 27ra^'"'[G(Mi + M2)]-'''^ (4) 

and assuming that the binary has a circular orbit one finds: 

/vise _ 37r« /Hdy" /fld\-^/^ , . 

/orb yrr^ \Rd) \ a ) ^ ' 

where q = M2/A/1 is the mass ratio, a is binary separation, 
/orb = 1/forb and /vise = lAvisc- From eq.0 it follows 
that for fixed viscosity parameter and relative disk thickness 
and size, the viscous time scale is directly proportional to 
the orbital frequency of the binary system, as intuitively 
expected. 

To estimate the expected range of values of /viae //orb 
one needs to know the disk size and thickness. The thickness 
of the standard Shakura-Sunyaev disk is 

I « 1.9 . 10- a-/^° 7??/- M-f L^f ° R\i' (6) 

where Mi = 1.4 Mi. 4 is the mass of the compact ob- 
ject in solar units, Ro — 10® -Riokm cm is defined so that 
Lx = GMiM/Ro, Lx = lO^'' L37 erg/s is X-ray lumi- 
nosity, and Rd — 10^" Rio cm is the outer disk radius. 
The disk thickness predicted by this relation varies from 



H/R ~ 0.02 for the most compact and low luminosity sys- 
tems to H/R ~ 0.03 — 0.04 for sources in the upper part 
of the Table Irradiation of the disk by the X-rays from 
the vicinity of the compact object, although does increase 
significantly the surface temperature of the outer disk, has 
a little effect on its mid-plane temperatu re, due to large op- 
tical depth of the Shakura-Sunyaev disk jLvutvi fc Sunvaievl 
Il976l) . Correspondingly, account for irradiation effects does 
not change the above estimates significantly. With plausi- 
ble values of two other parameters in eq.JSJ, 01 ~ 0.5 and 
Rd/a ~ 0.3 — 0.6, the standard theory predicts: 

(h^) - 0.005 -0.07 (7) 
4.1 Break frequency 

Power density spectra of all but one LMXBs from our sample 
have the low frequency break predicted in the simple qual- 
itative picture outlined in the Introduction. Fig. [7| shows 
the dependence of the best fit value of the break frequency 
on the orbital frequency of the binary system. There is an 
obvious positive correlation between these two frequencies 
broadly consistent with eq.Q. Remarkably, this correla- 
tion holds over 3 order of magnitudes in the orbital fre- 
quency, from compact binaries 4U 1820-303 and 1H1627- 
673 with Porb ~ 12 and ~ 40 min to extremely wide sys- 
tem GRS 1915+105 with Porb « 33.5 days. The value of 
/brcak//orb ranges from ~ 0.1 — 0.2 for long period binaries 
to ~ 1 — 2 for the more compact ones. The broad band spec- 
tra shown in Fig|^confirm that the breaks are indeed unique 
and well defined features in the power spectra. 

In the following we assume that the break frequency 
corresponds to the viscous time scale of the accretion disk, 
and they are related by a simple linear relation: 
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Figure 7. Relation between the PDS break frequency and the orbital frequency of the binary system. The shaded area shows the 
dependence /vise vs. farh expected for the Shakura-Sunyaev disk, assuming a = 0.5, /3 = /broak/Zvisc = li <? = 0.5 and Rd/a, = 0.4. 
Straight solid and dashed lines are predictions for larger values of the disk thickness H/R, as indicated by the numbers on the plot. The 
two points for Cir X-1 correspond to the observed orbital period of Porb ~ 16.6 d (open circle) and corrected for the orbital eccentricity 
assuming e = 0.7, as described in the section [4.2l ffiUed circle). 



^visc — /vise — /5 /break 



(8) 



Note that this simplified picture does not take into account 
complexity of the power density distribution near the viscous 
time scale. 

With this definition the break frequency is related to 
the orbital frequency via: 



/break _ ^ StTQ 

/orb VI + a 



4.2 Cir X-1 



(f)( 



-3/2 



(9) 



Cir X-1 appears to deviate from the general trend in Fig|2| 
having the break/viscous time scale by a factor ~ 10 shorter 
than expected given its orbital period. 

The significant orbital modulation of the X-ray activ- 
ity s uggests a highly ec c entric orbit in this b inary sys- 
tem iMurdin et al.lll980l) . Ijohnston et alj 1^9^ proposed 
a model, where the system consists of a neutron star on a 
highly eccentric, e ~ 0.7 — 0.9 orbit around a ~ 3 — 5 Mq 
sub-giant companion. In this model the donor star fills its 
Roche lobe and the mass transfer occurs only during the 
periastron passage. Outside the periastron the donor star is 
detached from its Roche lobe surface and the mass transfer 
in the binary stops. 

In such a system the disk radius would be defined, to 
first approximation, by the minimal separation between the 
components, amin = (1 — e) a. As the eq.l|SJ was derived for 
a circular orbit, the following substitution should be made 
in interpreting the Cir X-1 data: 



Porb -> Porb (1 - ef^"" 

Plotted in Fig|7]are two point for Cir X-1. The open cir- 
cle corresponds to the observed orbital period of the source, 
the filled circle corresponds to the orbital period corrected 
for the eccentricity of the orbit assuming e = 0.7. With this 
substitution the consistency with other sources is restored. 



4.3 Viscous time and disk thickness 

Although the break frequency does increase linearly with the 
binary orbital frequency, the ratio /broak//orb ~ 0.2 — 2 is 
notably larger than predicted for the Shakura-Sunyaev disk, 
( /vise/ /orb ~ 0.005 — 0.07. The latter range of values is 
shown as the hatched area in the Fig. |7| Obviously, larger 
values of /vise/ /orb iniply that the disk viscous time as traced 
by the position of the break on the PDS is by a factor of 
> 5 — 10 shorter than predicted by the theory. As it follows 
from eq.l|5J and @, there are several parameters affecting 
the disk viscous time, of which the strongest dependence is 
upon the disk size and thickness. 

The theoretical and observational constrains on the disk 
size are summarized in Fig|S| It has two well known the- 
oretical limits. Due to the angular momentum conserva- 
tion, the disk radius can not be smaller than the circular- 
ization radius. The r esults of the numerical calculations of 
iLubow fc ShrJ (1197511 can be approximated by 



Rciv 



0.074 



1 + g 



(11) 



a (1 



(10) 



This approximation is accurate to ~ 3% in the range of the 
mass ratios 0.03 < g < 10. The upper limit on the disk size is 
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Riven by the tidal truncation radius jPapaloizou fc Pringlj 
19771: llchikawa fc Osakilll994l) . In the case of small pressure 
and viscosity its value is close to the radius of the largest 
non- intersecting periodic orbit (iPaczvnskL1977.'l . The latter 
can be approximated by: 

0.270 , 0.239 

(12) 



= 0.112 



1 + g 

3% : 



accurate to ~ 3% in the range of 0.06 < 9 < 
10. These two limits on the disk radius are shown 
as thick dash-dotted lines in Fig|H| On the observa- 
tional side, there is a number of spectroscopic and 
phot ometric determin ations of the dis k radius in CVs 
ilHessman 1988; Hcssmaii fc Hoop 1990: Rutten et alj|l992t 



iHarlaftis fc MarshlEgg'dlHarrop-Alliri fc WarneJll99gi and 



fewer in LMXBs ilOrosz fc Kuulkerd Il999l: IShahbaz etlll 
|2004| : iTorres et ai]|2004l : IZurita et aljbooof) . These data in- 



dicate (Fig|5J that for a steady disk the disk radius is close t o 
the tidal truncation radius as defined by Paczvnskj il977tl . 
with possible exception of the small-q systems f section 14.411 . 

With these values of the disk size and accepting a ~ 
0.5 as a plausible value of the viscosity parameter, the disk 
thickness of 



H/R>Q.l a-'^''l3 



1/2^-1/2 



(13) 



is required in order explain observed values of the break 
frequency. 

This conclusion relies on the association of the break 
in the power density spectra with the viscous time scale of 
the disk and as long as the assumption /3 = /brcak//visc ~ 1 
is valid, it is very robust. Indeed, in order to describe the 
data with H/R ~ few x lO"'^, as predicted by the stan- 
dard theory, one would need to increase dramatically the a- 
parameter, up to implausibly high values of a ~ 5 — 800. Al- 
ternatively, short viscous time of the disk could be achieved 
by reducing the Roche-lobe filling factor of the disk, down 
to Rd/a ~ 0.005 — 0.05. This is significantly smaller than 
the circularization radius and contradicts to both theoreti- 
cal expectations and observations of the disks in CVs and 
LMXBs (FiglHl. 

The latter possibility can be also interpreted that only 
the very inner part of the accretion disk, R/a ~ 0.005 — 
0.05, contributes to the observed variability of LMXBs in 
X-rays. This can not be excluded a priori. We note however, 
that in this case the sharpness of the breaks observed in the 
power spectra of at least some sources (e.g. Sco X-1, EXO 
0748-676, etc) implies that the transition from the (outer) 
region of negligibly small amplitude of the M perturbation 
to the (inner) region of relatively large ones should occur in 
a rather narrow range of radii. 

4.4 Viscous time in wide and compact systems 

As is evident from Fig|7] there is a significant dichotomy in 
the value of /break //orb ratio between long- and short-period 
binaries. This is further illustrated by Fig|^ where the ratio 
/break //orb is plotted agaiust the binary orbital period and 
the mass ratio. The average values of /break//orb are « 0.23 
and ~ 1.7 for the wide and compact systems correspond- 
ingly. The sense of the difference is that the compact sys- 
tems have systematically shorter, by a factor of ~ 8, viscous 
time expressed in the orbital periods than the wide ones. 
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Figure 8. Theoretical and observational constrains on the disk 
radius. The solid curve shows the effective radius of the Roche 
lobe of the primary. Two dash-dotted curves show the circular- 
ization and tidal truncation radius. The dashed line is the radius 
of 3:1 tidal resonance. The open circles and triangles are measure- 
ments and lower limits for CVs in the outburst state. The small 
filled circles show the disk radius in LMXBs. The reference for 
these optical data are given in the text. 



the boundary lying at Porb ~ 10 hours or g ~ 0.3. The fact 
that this bimodality occurs between compact and wide sys- 
tems or, equivalently, (with the exception of GRS1915-f 106) 
between small- and large-q binaries suggests that it may be 
caused by the excitation of tidal resonances in the accretion 
disk. 

The phenomenon of tidal resonances in ac cretion disks 
is we l l studied in the context o f dwarf novae iWhitehursti 
119881: IWhitehurst fc Kind [l99ll: iLubowl Il99l . The reso- 
nance occurs if the angular frequency of the orbital motion 
of the particle in the disk is commensurate with the angu- 
lar frequency of the orbital motion of the secondary. For a 
Keplerian disk, the location of the resonant orbits is given 
by 

^^(A)^^%l + .)-/3 (14) 
a \mj 

where k, m are integers (e.g. IWhitehurst fc King|ll99j) . An 
obvious condition for a resonance being excited is that 
this radius lies within the accretion disk. From analyti- 
cal studies and numerical simulations it is known that the 
strongest resonance occurs at t he lowest order commensu - 
rability, : i^orh = 2:1 (e.g. IWhitehurst fc Kin'3ll99ir) . 
However the 2:1 resonant radius lies within the accretion 
disk only at extreme values of the mass ratio, g < 10~^. The 
next strongest resonance, most important in the context of 
binaries with not too extreme values of the mass ratio lies at 
the 3:1 commensurability. The higher order resonances are 
probably not excited in the accretion disks. The dependence 
of the 3:1 resonance radius on the mass ratio q is shown by 
the dashed line in Fig|H| confirming the well known fact that 
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Figure 9. Dependence of the ratio /brcak//orb on the orbital period of the binary system {left) and its mass ratio (right). 



the 3 : 1 resonance can be excited in the systems with the 
mass ration g<0.3 — 0.4. Note that the precise value of the 
threshold q depends on the definition of the tidal truncation 
radius, due to small angle between the two curves in Fig|H| 
Under the assumption that i?tidai = 0.9Rl it is g ~ 0.33. 
With the definition of i?tidai used in Fig|H|the threshold value 
is slightly larger. 

Numerical simulations of the non-linear stage of the in- 
stability have shown that excitation of tidal resonances re- 
sults in a significant asymmetry of the accretion disk and 
causes its precess i on in the reference frame of the binar y 
l|Whitehurstill988l: IWhitehurst fc KindlToQll: lLubowlll99ll) . 
Additionally, the accretion disk is truncated near the res- 
onant radius. This explains the values of the disk radii 
smaller than the tidal truncation radius observed in CVs 
and LMXBs with small mass ratios (FiglUl. 

To conclude, with exception of GRS1915-I-105, dis- 
cussed in section 15. 2. II the excitation of the 3:1 inner Lind- 
blad resonance provides a natural explanation of the di- 
chotomy in the viscous time scale between wide and compact 
systems. 



5 DISCUSSION. 

5.1 Thickness of the disk 

5.1.1 Semi-thick disk or coronal flow? 

The values of the viscous time scale of the disk inferred 
by the position of the low frequency break in the power 
spectra of LMXBs require rather large disk thickness, 
H/R> 0.1 a^g/^/J-^/^ by a factor of >3 - 5 exceed- 
ing the prediction of the standard theory of the accretion 
disks. Taken at the face value, H/R> 0.1 appears to agree 
with the statistics of the eclipsing systems among LMXBs 



fe.g. lMilgromlllQTStlNaravan fc McClintockl20oi) . However, 
large values ol H/R result in small optical depth of the outer 
disk, which is inconsistent with numerous observations of the 
optical emission from accretion disks, as discussed below. 

We consider the structure of an accretion disk with 
some value of the vertical scale-height H/R ~ 0.1, with- 
out specifying the source of the additional energy dissipa- 
tion. Under the assumption of vertical hydrostatic equilib- 
rium and stationarity, both the mid-plane disk temperature 
and the density are steep functions of the disk thickness: 
r/Tvir oc {H/Rf, p oc {H/R)-^. Due to strong dependence 
of the Rosseland mean opacity on the density and tem- 
perature, typically as kn oc y»-3.5 ^j^g parameters 
range of interest, the optical depth of the disk tr = kapH 
is also steeply decreasing function of the disk thickness, 
tr oc {H/R)~"' with the power law index 7 changing from 
~ 12 to ~ 3 as the H/R increases (Fig llUl . Therefore, 
even a moderate increase of the disk thickness would lead 
to dramatic decrease of its optical depth. As a result, for 
H/R ~ 0.1, the accretion disk in the majority of LMXBs 
from our sample would be optically thin. This would con- 
tradict to a number of observational facts, mostly from the 
optical band, proving existence of the optically thick outer 
disk in LMXBs. Indeed, the free- free optical depth in the 
optical V-band equals 



f"//(V-band) 



(15) 



3.8 • 10" 



' In f -^—] Mfr M;1'^ a^l ( -) R 

V5500 kJ \ rJ 0.1 



-1/2 



where {H/R)o,i is the disk geometrical thickness normalized 
to 0.1. The above approximation is accurate to ~ 10% in the 
10* < T < 10* K temperature range, the logarithmic term 
due to the Gaunt factor varies from 2.9 to 7.5 for the temper- 
ature in the range 10^ — lO'^ K. For the parameter of binaries 
from our sample the optical depth in the V-band would be 



Viscous time scale and Lindblad resonances in LMXBs 13 



The vertical and radial column density of the coronal 
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Figure 10. Dependence of optical depth of the disk at _R = 5-10^" 
cm on the H/R ratio for M = 5 ■ 10^'' g/s, M = 1.4 M©, a = 
0.5. The solid line shows the opt ical depth based on the OPAL 
opacities ilglesias fc RogersllQOd) , long dashed line - Kramer law 
(kji = 5T0^* pT~'^/'^ cm^/g), short dashed — free-free absorption 
in the optical V-band, dash-dotted - Thompson optical depth. 



Tff <^ 1, with the exception of sources with the largest M 
(Sco X-1, GX349-h2 and 4U1820) where t// ~ 1 - 10. The 
low optical depth in the V-band would conflict with numer- 
ous spectroscopic observations of the optical emission from 
the LMXBs, suggesting existence of the continuum emission 
originating f rom the o ptically thick o uter accretion disk (e.g. 
iHvnes et al]l2 002: Cha tv et al1l2003ll . On this basis the pos- 
sibility of the semi-thick, H/R ~ 0.1 outer disk in LMXBs 
can be excluded. 

Instead, one can envisage a two phase accretion flow 
with a Shakura-Sunyaev-like disk in the mid-plane sur- 
rounded by the tenuous optically thin coronal flow with 
the aspect ratio H/R ~ 0.1 — 0.2 and temperature of 
T ~ 10~^Tvir. In order to explain the prominence of the 
low frequency break in the power spectrum, a significant 
fraction of the accretion has to take place in this corona 

Mcorona ~ Mdisk (16) 

The temperature and density in the coronal flow are: 



Tr^ 1.4- lO"" Mi.4 (—) /in' [K] 
V it / 0.1 



3.6 • 10" /e Mi7 M-y^ an 



°'\r)o.i^'' 



(17) 



(18) 



where fc = Mcorona/M. As the surface mass density depends 
on the thickness of the disk as E oc (H/R)-'^, such a corona 
will contain a small fraction of the mass of the accreting 
matter: 

'<10"' (19) 



fc 



assuming /c ~ 0.5. 



flow are \cm 



nH^3.6- 10^'' fc A/i7 M'^^ a, 



ni? « 3.6 



10'''' fc Mi7 M-y^ a, 



°'K:r)o.i 



1/2 



-^11 



1/2 



(20) 



The plausible mechanism of formation of the coronal 
flow is the disk evaporation proc ess originally proposed by 
iMever fc Mever-Hofmist'eil (^993) for dwarf novae. In X-ray 
binaries the physics of the disk evaporation is significantly 
affected by the illumination and Compton heating and cool- 
ing of the disk-corona system by the X-ray emission pro- 
duced in the vicinity of the relativistic object. Although the 
the self-consistent treatment of the problem is yet to be 
done, the likely net effect of the X-ray irradiation is to in- 
crease the fraction of M in the coronal flow ( Me ver et alJ 
I2OOOI : Ijimenez-Garate et alJl2002L F.Meyer, private commu- 
nication). 

5.1.2 Other evidence of the coronal flow 

An independent evidence of the existence of a diffuse ionized 
gas above the accretion disk plane is provided by observa- 
tions of LMXBs with high inclination - the ADC sources and 
dippers. Bas ed on observations of partial X-ray eclipses in 
4U1822-371, IWhite et all dlQSlT) concluded that the com- 
pact X-ray source in this system is diffused by a large 
moderately Compton thick highly ionized corona located 
above and below the accretion disk. Detailed modeling of 
the eclipse li ght curves in the AD C sources 4U1822-371 and 
4U2129+47 iwhite fc Hoi^llQSall yielded estimates of the 
radial extent of the accretion disk corona in these sources, 
~ (1 - 2) • 10" cm and i?c ~ (3 - 6) • 10^°. These num- 
bers are comparable with the value of the tidal radius in 
the systems, J?tidai ~ 3.9 • 10^" cm. The partial nature of 
the eclipses in these sources, with the eclipsed flux at the 
level of ~ 50% of the uneclipsed value indicated that the 
corona had non-negligible radial optical depth, tt ~ 1, in 
agreement with the estimate of ea. ()2()^ . 

High resolution spectroscopic observations of LMXBs 
with Chandra and XMM-Newton gratings are reveal- 
ing com plex absorptiony emission features in their X-ray 
spectra jC ottam et al.' 2001; "jimene z-Garate et alJ l2003l : 
iKallmanet aL.2003. : .Boirin et aL.2004h . These features are 
mostly pronounced in the systems with high binary inclina- 
tion angle and suggest the presence of tenuous photo-ionized 
plasma orbiting the compact object above the orbital plane 
of the binary. 

From the analysis of the recombination emission lines 
of H- and He-like ioii s of Q , Ne and Mg in EXO0748-676 
I.Iimenez-Garate et al.l (|200 .3h constrained the density and ra- 
dial extend of the photoionized plasma, Nc ~ 10^^ — lO'^^ 
cm"^, 7-10^ <i?<6-10^ cm and the temperature T < 2 ■ lO'^ 
K. This range of the parameters is consistent with the num- 
bers computed from the above formulae for EXO0748-676 
and assuming the aspect ratio of the coronal flow H/R = 0.1: 
the tidal truncation radius iitidai ~ 4.8 • 10^" cm, density in 
the coronal flow A'' ~ 2.7 ■ 10'^''/c cm^'', and temperature of 
T~ 3 - lO'^ K. 

Based on the observations of the narrow absorption lines 
of the H- and He-like Fe in the spectrum of 4U1916-053, 
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iBoirin et alJ ll2004l) estimated the ionization parameter and 
the column density of the absorbing gas: log(^) « 3.9, NH ~ 
2 ■ lO'^^ cm"'^. From these numbers one can estimate the 
radial extent and the density of the absorbing gas, R ~ 
1.6 • 10^" cm, iV ~ 1.3 ■ 10^^ cm~^. The tidal radius in this 
system is « 2 • 10^° cm, in agreement with the value derived 
from the XMM-Newton data. The density in the coronal 
flow with H/R = 0.1 is AT ~ 6.7 ■ lO^Vc cm'^ The latter 
number is somewhat larger than derived from the XMM- 
Newton data. This can be explained by smaller inclination 
angle in this source resulting in lower density along the line 
of sight. Indeed, unlike EXO0748-676 showing deep X-ray 
eclipses, 4U1916-053 shows only X-ray dips. 

5.1.3 M variations due to geometrically thin disk 

In this picture, the M variations in the geometrically thin 
Shakura-Sunyaev disk should give rise to the second very 
low frequency power law component, revealing itself below 
the break frequency. Such power law component is indeed 
observed in the power spectra of some of the sources. Most 
obvious its presence is in Cir X-1 and 4U1916-053 and, to 
less extent, in Sco X-1 (FigEJ. Remarkably, the slope of the 
power spectrum before and after the break are very similar, 
suggesting similar nature of the processes causing the vari- 
ability. Furthermore, the second break should be expected 
at the low frequencies, corresponding to the (larger) vis- 
cous time of the geometrically thin disk. The frequency of 
the second break is described by eq.Q with the H/R given 
by eq.JSJ. An evidence of the second break may be seen 
in the power spectra of Sco X-1 and, possibly, of 4U1916- 
053. In the former case the break frequency equals fy,.^i,^y_ = 
■ 10"* Hz and the ratio fiH^Jf orb = 5^1 ■ 10"^ This 
value is within the range predicted for standard Shakura- 
Sunyaev disk with the aspect ratio of H/R ~ few • 10~^, 
eq.Q. The value of the break frequency in the case of 
4U1916-053 is probably /^^^^^ ~ lO"'' Hz, ahhough this 
number is less secure than that for Sco X-1. If correct, this 
is significantly, by ~ 2 — 3 dex smaller that expected in the 
above scenario. 



In the conclusion of this section we note that due to 
power law dependence of the viscous time scale on the disk 
radius, tvisc oc r^^^, the radial drift velocity of the matter 
in the accretion flow should be increased throughout the 
most (~ 80 — 90%) of the radial extent of the accretion disk. 
Therefore the bulge at the outer edge of the disk, result- 
ing fr om the disk-stream interaction (e.g. Hollicr & Mason 
119891) is insufficient as it is located at the outer edge of the 
disk. For the same reason this analysis based on the viscous 
time of the disk probes the outer disk and is insensitive to 
the conditions in its innermost parts, e.g. i?< 10* -Rg, which 
contribution to the total viscous time is negligibly small. 

5.2 Tidal resonances and viscous time scale of the 
disk 

The dichotomy in the viscous time between long- and short- 
period systems is naturally explained by excitation of the 3:1 
inner Lindblad resonance. This suggestion is supported by 



detection o f superhumps in several LMXBs with sm all mass 
ratios (e.g. iHaswell et aljEoOll : IZurita et al.ll2002l) . Impor- 
tantly, the superhumps are detected not only in transient 
sources but also in one or two persistent L MXBs, including 
one s ource from our sample, 4U1916-053 dCallanan et alJ 

The specific mechanism by which the tidal instability is 
affecting the viscous time of the disk is, however, unclear. 
It is also unclear if the coronal flow in the small-q system 
has indeed H/R ~ 0.3 and T ~ O.lTvir as it is formally 
required by the observed fvisc ~ Poih. Several possibilities 
can be mentioned: 

(i) mass transfer in tidal waves excited in the disk by the 
tidal forces. It has been shown theoretically and confirmed in 
numerical simulations, that in the non-linear r egime strong 
tidal waves are excited in the accret ion disk dSpruit et alJ 
119871: lLubo'^ll99ll: iTruss et al.ll2002^ . These waves, propa- 
gating in the disk with the speed significantly exceeding the 
radial drift velocity of the matter, can decrease the effective 
viscous time scale. If correct, the aspect ratio and coronal 
temperature in the small-q systems may be not significantly 
different from those in the large-q ones. 

(ii) hea ting of the outer disk due to resonant tida l in- 
teraction JWhitehurst fc Kindfllill : iTruss et alj|2002l) can 
increase the temperature of the disk and/or coronal ffow 
and, correspondingly, decrease the viscous time scale. As 
was mentioned above, in order to explain the observed vis- 
cous time scale the coronal flow temperature of T ~ 0.1 Tvir 
is required. 

(iii) non-trivial definition of the viscous time for an eccen- 
tric and precessing disk. In the fully developed instability the 
particles in the disk move on non-circular trajectories with 
large eccentricity. Therefore, the average velocity of the ra- 
dial drift of the disk particles can be significantly different 
from the value predicted by the standard disk theory for a 
circular Keplerian disk. 

We note, that the truncation of the accretion disk at the 
3:1 resonance radius (FiglSJ can be excluded as the cause of 
the shorter viscous time scale in the small-g systems, be- 
cause the change of the outer disk radius by a < 20 — 30% is 
insufficient to explain the observed decrease of the viscous 
time by the factor of ~ 10. 

5.2.1 Case of GRS1915+105 

GRS1915-(-105 clearly stands out among other small-q sys- 
tems (Figl^ right pan el). Although it has small mass ra- 
tio, q = 0.058 ± 0.033 jHarlaftis fc Greineilliooi) . its ratio 
/brcak//orb IS close to the valuc fouud in the binaries with 
g>0.3. The most plausible explanation of such behavior is 
related to the transient nature of the source. We outline two 
possibilities: 

(i) the tidal instability is known to have a rather long 
gro wth time, of the o rder of > 10^ binary orbital periods 
fe.g. lWhitehursiigSgH . For GR^1915-H05 Porb ~ 33.5 days 
and 10^ Porb corresponds to ~ 10 years. Therefore the time 
passed since the onset of the outburst (~ 10 years) was 
insufficient for the tidal instability to fully develop. In this 
scenario one would expect an increase of the break frequency 
in the course of the outburst, unless the instability growth 
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time is significantly longer than ~ 10 years. However, no 
statistically significant difference in the break frequency has 
been found between the first and second halves of the data. 
This fact speaks against the above suggestion. 

(ii) due to the large size of the system (~ 8 ■ 10^^ cm), 
the irradiation effects are insufficient to ionize the outer 
disk in GRS1915+105 which is in the cold, low-viscosity 
state (H.Ritter, private communication). Therefore its outer 
boundary is located near the circularization radius rather 
than has expanded to the tidal radius, as is the case for 
the persistent sources. Therefore there is no accretion disk 
at the 3:1 resonance radius (Fig|HJ and the tidal insta- 
bility is not excited. In this respect the present outburst 
of GRS1915-I-105 might be similar to normal outbursts in 
dwarf novae. As well known, for the same reason the su- 
perhumps in dwarf novae are observed (tidal instability is 
excited) only in during th e superoutbursts but not during 
the normal outbursts (e.g. lOsaki fc Meveill20o3) . 



5.3 Broad low frequency QPO features 

For several sources broad QPO-like features are obvious 
in the power density spectra. Most apparent these fea- 
tures are in the case of Cyg X-2 and GRS1915-I-105 
(Figl^l I3J . Among plausible mechanisms of appearance of 
quasi periodic variability in X-ray binaries on these time 
scales are the radiation-d riven warping and p recession of 
the accretion disk (e.g. iMalonev. Beeelman fc Prineie 1996) 
and mass-flow oscillations cause d by diffusion instability 
jMever fc Mever-Hofmisteilll99d) . Both types of oscillations 
are caused by the irradiation of the accretion disk by the 
primary and their characteristic time scale is of the order 
of the viscous time scale of the disk. Accordingly, these fea- 
tures should appear near the break frequency in the power 
density spectrum. 

Another possible mechanism of appearance of broad low 
frequency QP Os near the viscous frequency of the di sk was 
suggested bv IVikhlinin. Churazov fc Gilfanovl il994) . This 
mechanism is related to the anti-correlation of the M per- 
turbations in the accretion flow on the time scales t ~ tvisc. 
The ant i- correlation appears in the presence of the steady 
mass supply through the outer boundary of the accretion 
flow and is a direct consequence of the energy conservation 
law. In this respect the broad QPO features observed near 
the viscous time scale of the outer disk might be similar to 
the broad QPOs detected near the high frequency break of 
power density spectrum of C yg X-1 in the hard s tate, at 
frequency of ~ 0.04 - 0.07 Hz ijvikhlinin et alJll994ll . In the 
latter case the values of the break frequency and of the QPO 
centroid frequency are defined by the viscous time scale of 
the inner hot fiow, R < 10^ Rg , responsible for the produc- 
tion of the hard Comptonized emission in the low spectral 
state of black hole candidates. 



(i) As the viscous time is the longest time scale of the 
accretion disk, the M variations produced in the disk due 
to various instabilities should become independent of each 
other on the time scales longer than the viscous time, 
i.e. (^M{t)M{t + t)^ — + at r>tvisc Correspondingly, the 
power density of M perturbations produced in the disk and 
of the observed X-ray fiux variations should become inde- 
pendent of the frequency at /^tvisc- 

(u) Using archival data of RXTE/ASM and EX- 
OS AT/ME we studied X-ray variability of persistent 
LMXBs in the / ~ 10"* - 10"^ Hz frequency range. In 
power density spectra of 11 sources out of 12 satisfying our 
selection criteria, we found a very low frequency break. The 
spectra approximately follow oc power law above 

the break and are nearly fiat below the break (Figs. EHSJl. 
In some cases (Sco X-1, Cir X-1 and 4U1916-053) a sec- 
ond power law component with the same slope appears at 

/ ^ .fbrcak. 

(iii) In a very broad range of binary periods, from Porb ~ 
11 min in ultracompact binary 4U1820-303 to Porb ~ 33.5 
days in GRS1915-f 105, the break frequency correlates with 
the orbital frequency /orb = l/Porb of the binary system 
(Fig|7|l, in good agreement with the theoretical prediction 
for the viscous frequency /vise = t^isc (sfllSl 

(iv) Assuming that the low frequency break is associated 
with the viscous time scale of the disk, /break ~ l/^visc, the 
measured values of the break frequency imply that the vis- 
cous time of the disk is by a factor of > 10 shorter than 
predicted by the standard theory f section 14.31 Fig|7J. Taken 
at the face value, this requires the relative height of the 
outer disk H / R>0.1. Motivated by the low vertical optical 
depth of such semi-thick accretion fiow we propose instead, 
that significant fraction of the accretion occurs through the 
coronal flow above the standard geometrically thin Shakura- 
Sunyaev disk (section l5.1.1l Fie llOH . The aspect ratio of the 
coronal flow, ///P > 0.1, corresponds to the gas temperature 
of r>O.Oirvir. The corona has moderate optical depth in 
the radial direction, tt ~ 1, and contains < 10% of the total 
mass of the accreting matter. The estimates of the tempera- 
ture and density of the corona are in quantitative agreement 
with the parameters inferred by the X-ray spectroscopic ob- 
servations by Chandra and XMM-Newton of complex ab- 
sorption/emission features in the LMXBs with large incli- 
nation angle. 

(v) We found a clear dichotomy in the /broak/./orb between 
wide and compact system, the accretion flow in the compact 
systems having ~ 10 times shorter viscous time expressed 
in the orbital periods of the binary, than in the wide ones 
(Fig|UJ. The jump in the ivisc/Porb ratio occurs at the binary 
system mass ratio of q ~ 0.3. This strongly suggests that the 
dichotomy between small-q and large-q systems is caused by 
the excitation of the 3:1 Lindblad resonance in the small-q 
systems (section [4.411 . 



6 SUMMARY 

We studied the low frequency variability of low mass X-ray 
binaries in order to search for the signatures in their power 
density spectra related to the viscous time scale of the disk. 
Our results can be summarized as follows: 
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APPENDIX A: POWER SPECTRA FROM THE 
THE ASM DATA 

To calculate power density spectra from ASM light curves 
we used the method based on the auto-correlation function 
(cf. lEdelson & Krolik 198^ . 

We begin with consideration of the equally spaced data. 
The measured time series is Xk and equals the number of 
counts detected in the fc-th time bin corresponding to the 
time interval of t = (fc— 1) At-^fcAt, where index k runs from 
to A'^ — 1, N is assumed to be even. The discrete Fourier 
transform of the time series is defined as : 



Xk exp 



(Al) 



■^j — / , 

fe=0 

where j changes in the range 1 N 12. The power density 
in the ji-th frequency bin, which center is = j/T Hz, 
expressed in units of rms^/Hz (rms is fractional rms of vari- 
ability) equals: 

p ^ 2|i£^ ^ 2!i,p 



iVph (r) (r)'r 



(A2) 



where TVph = ^^q^ is the total number of counts in the 
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time series, T — NAt is its total duration and (r) — Np^/T 
is the average count rate. 

Combining eas lAll and I A 21 one finds 



JV-l 



{rfT 



AoN + 2'^{N- 



- 1) (Ai) cos 



iV 



(A3) 



where (Ai) is discrete auto-correlation function: 

(Ai) = {xk Xk+i) (A4) 

The eg I A3I expresses the well known fact that power equals 
the cosine transform of the auto-correlation function. The 
noise level in the power density spectrum equals: 



EJV-l 
fc=0 



(A5) 



where Ufc is the error of Xk- For Poisson errors in x^, cr^ — Xk 
and Pnoise = 2/ (r), in a greement with the conventional form 
(e.g. iLeahv et al] ll98^ . 

The ASM time series consists of the count rate mea- 
surements r-k, [cnts/sec], at unevenly spaced moments tk, 
fc = 1, A'asm- The values of rk are averaged over the time 
interval (~ 90 sec) much shorter than the distance between 
adjacent bins (~ 90 min). As before, the auto-correlation 
function Ai is defined on the grid of A'' time bins, with the k- 
th bin corresponding to the time interval t = (fc— l)At-^fcAt, 
where index k runs from to A — 1. The {Ai) is computed 



{A{) = {{n-{r))(r,-{r))) At' 



(A6) 



where the averaging is performed over all pairs i > j, 

for which Atij = ti — tj falls in the /-th time bin of the grid. 

The time interval At and the duration of the time series 
T = NAt are parameters of the transformation and define 
the frequency range covered by the output power density 
spectrum. Their choice depends on the effective Nyquist fre- 
quency of the time series and its duration. Appropriate for 
the parameters of the ASM light-curves are the following 
values: At ~ 5 • 10^ 4- lO" sec and T ~ 10^ 10** sec. Ad- 
vantage of this method is that the correctness of the chosen 
values of At and T — NAt can be easily checked using the 
number of pairs used to compute the autocorrelation 

value at the given time scale I At and the uniformity of these 
numbers across the full range of the time scales, from At to 
T = NAt. 

With the auto-correlation computed according to 
eg IA6I the power density spectrum can be easily computed 
using eg IA3I As usual in the time series analysis, the en- 
tire time series is divided into n segments and the power 
spectrum is computed for each segment separately. Their 
average and dispersion give estimates of the average power 
density spectrum and its uncertainty (under the assumption 
of stationarity) . 

If the errors for are known, the noise level in the 
output power spectrum can be computed as 



P 

no 



(A7) 



(r) Aasm 

where ak is the error of rk (in units of cnts/sec) and Aasm 
is the number of ASM measurements used to compute one 
power spectrum. In practice, the fiux measurement errors 
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Figure AI. Results of simulations. The power density spectra 
of ASM light-curves simulated as described in the text. The solid 
lines show the model power density spectra used to generate faked 
ASM light curves: power law = 10~*/~^'^ with low frequency 
break at 10"^ and 10"'' Hz. The model with /break = 10"'' Hz 
and corresponding power density spectrum are shifted upwards 
by the factor of 4 for clarity. 



given in the ASM light-curves are slightly underestimated, 
therefore the noise level computed with eg I A 71 is somewhat 
imprecise (Grirmn ct_al. 2002). 

The performance of the method described above was 
verified in simulations. The initial light curve with the time 
resolution of ~ 100 sec was computed assuming the power 
law spectrum with slope of 1.3 and low frequency break at 
/ — 10~® Hz and / — lO"'^ Hz and random Fourier phases. 
From this light curve the ASM time series was simulated us- 
ing the measurements times tfe of the real ASM light curve of 
Cyg X-2. The obtained light curve was randomized assum- 
ing Gaussian errors with the standard deviation equal the 
real ASM measurement error for the given measurement rk- 
The obtained time series was analyzed with the same code as 
used for the analysis of the real ASM data with parameters 
At = 10* sec and T — 5 ■ lO' sec. The noise level estimated 
using eg |A7l was subtracted from the power density spectra. 
The results of simulations are shown in Fig lAll along with 
the model power density spectra used to generate the initial 
high resolution light curves. 



